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Abstract 

The influence of tlie impurity on the resonant Andreev reflection through a normal-metal/carbon- 
nanotube/superconductor system is studied theoretically. It is found that the resonant Andreev 
reflection depends on the strength of the impurity and the length of the armchair nanotube. The 
impurity which breaks the electron-hole symmetry of the nanotubes greatly reduces the resonant 
Andreev reflection. The symmetry broken depends distinctly on the impurity strength. The 
impurity effects on the Andreev reflection current at different bias are also studied. 
PACS number(s): 72.80.Rj, 73.23. Ad, 73.63.Fg 



Carbon nanotubes have been the subject of an increasing number of experimental and the- 
oretical studies due to their quasi-one-dimensional structure and unique electronic property.— 
The perfect carbon nanotube is predicted to be either metallic or semiconducting sensitively 
depending on its diameter and chirality, which is uniquely determined by the chiral vector 
(n, m), where n and m are integers.i^iM Experimental and theoretical studies have indicated 
that the electronic and transport properties of carbon nanotubes can be substantially modi- 
fied by point defects such as the substitutional impurities.— >^»2»^»2ii2iii One of recent interests 
concentrates on the electron transport through hybrid nanotube system. Experiments about 
some hybrid systems including nanotube-based magnetic tunnel junctionsi^ and supercon- 
ducting junctionsi^ii^ have been successfully fabricated. The electrical transport transport 
about the carbon nanotube quantum dot in the Kondo regime coupled to a normal and a su- 
perconductor has also been reported. ^^ The theoretical investigation of transport properties 
of these hybrid nanotube devices is of great importance, not only for their basic scientific 
interest, but also aiming at the design of novel nanodevices. 

The resonant Andreev reflections in the superconductor/carbon-nanotube devices 
has been theoretically studied.— The proximity effect in superconductor/carbon- 
nanotube/superconductor (S/CNT/S) tunnel junctions has also been studied theoretically.— 
However, the effects of the impurity, such as the substitutional boron (nitrogen) and the va- 
cancy, in these systems are not considered. In this paper, the effects of the substitutional 
boron (nitrogen) and the vacancy on the resonant Andreev reflection in the hybrid normal- 
metal/carbon-nanotube/superconductor (N/CNT/S) system are theoretically studied. In 
such a system, the specific molecular orbital plays an important role. By combing standard 
nonequilibrium Green's function (NGF) techniquesti^*^^^ with a tight-binding model ,^^^^^ 
we have analyzed quantum transport properties of the N/CNT/S system with the substi- 
tutional boron (nitrogen) and the vacancy. The Andreev reflection through the finite-sized 
carbon nanotube depends on the impurity strength and the tube length. The substitutional 
impurity, such as the boron and nitrogen, reduces the resonant Andreev reflection greatly. 
However, the impurity with very large strength, such as the vacancy, does not reduce the 
resonant Andreev reflection. The dependence of the Andreev reflection current with the 
gate voltage is also studied. 

We assume that the system N/CNT/S under consideration is described by the following 



Hainiltonian: 

H = Hl + Hr + Htube + Ht, (1) 

where 

Hl = ^{^tk - evL)a{j^^aL,ka, (2) 

k,cr 

Hr = 2_^ '^R,pO'R,paO'R,pcT + 2^[^*0/?,pi«/?,-pT + ^«iJ -pT^^'PJ-]' 
p,cr p 

Htube = ^ bo4aCja + H.C.] - ^ CVgcl^Ci^ + t/cJcQ, 
{i,j),a i,a- 

Ht = J2^t,al,^c,^ + H.c] + 5^[ti?e-^«-alj_^^c,. + H.c], 

k,a,i P>f,i 

where Hl describes the noninteracting electrons in the left normal- metal lead, ci\ f^^{aL,ka) 
are the creation (annihilation) operators of the electron in the left lead, and vl is the voltage 
of the left lead. Hr describe the right superconducting lead with the energy gap A. The 
nanotube Hamiltonian Htube is described by the tight-binding model with one vr-electron per 
atom at i site. The sum in i, j is restricted to nearest-neighbor atoms, and the bond potential 
7o = — 2.75e\^, which is used as the energy unit. This model is known to give a reasonable, 
qualitative description of the electronic and transport properties of carbon nanotubes .^^i^^ 
For simplicity, the intra-tube electron-electron Coulomb interaction has been neglected. Vg is 
the gate voltage which controls the energy levels in the CNT. The pointlike defect is defined 
by setting site energy equal to U at one of the sites of the unit cell, and various strengths 
represent typical substitutional impurities or vacancy. ^°. According to former tight-binding 
and ah inito calculations , ^°i^^ we set the strength f/ = 3, —5, and 10000 to simulate the 
substitutional boron, nitrogen and vacancy, respectively. H^ denotes the tunneling part of 
the Hamiltonian, and ti^n are the hopping matrix. It is convenient to introduce the 2x2 
Nambu representation in which the Green's function can be expressed by 

, V- AciT(r),4,(r')) (Cit(r),c,-|(r'))\ 

The retarded Green's function of the nanotube is calculated directly using a tight-bingding 
model via 
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g^{e) = I ^-^-'"=+*"" ^ I . (4) 



Using the standard NGF techniqu o^^i^^i^° , the Green's functions are obtained as 

g:,(.) = |fe(^))-' + 'fL + 5;^==r« + jf^^f (')]-', (5) 

>!'(<) = lfeW)-' + ^r, + ^-JlL=r„i-.. 

where F^, ^ are the appropriate hnewidth functions describing the couphng of the CNT to 
the respective leads. Here Gu and G12 are the retarded Green's functions of the CNT, which 
include the proper self-energy of the leads .^^i^^i^" Then the current and probability of the 
Andreev reflection are given by 

/a = ^ y delhie + cvl) - hie - et;z.)]T^(e), (6) 



T^{e)=Tr[r,G[,{e)r,Gll{e)], (7) 

where Jl^r denote the Fermi functions of the left and right leads, respectively. Clearly, the 
conventional tunneling is completely forbidden for V^ < A, and only the Andreev reflection 
exists. In the following numerical calculations, we discuss in detail the Andreev reflection at 
zero temperature in the case of V^ < A. We set (1) the temperature T = 0, (2) the voltage 
of the right lead vr = due to the gauge invariance, and carry out all calculations with 
A = 1, and F^ = F/j = 0.02 in units oi h = e = 1. 

In the following, the probability T^ and current I a of Andreev reflection for the N/CNT/S 
hybrid system are calculated. In Fig. 1, T4 is plotted as a function of the incident electron 
energy for nanotubes with different lengths. For the comparison, the conductance G of 
the N/CNT/N system with L = 3 is plotted in Fig. 1(a). The conductance peaks reflect 
the band structures of the finite-sized nanotubes, because the resonant states are close to 
the eigenvalues for small coupling F. In the vr-electron tight-binding model, the defect-free 
nanotubes have complete electron-hole symmetry with their Fermi levels at zero.— Then the 
resonant states are symmetric around the Fermi energy Ep = 0. The Andreev reflection 
probability for the N/CNT/S system with L = 3 are plotted in Fig. 1(b). The resonant 
Andreev reflection reflects the distribution of the resonant states of the nanotube systems. 
Although the resonant peaks are similar to those in a normal N/CNT/N system, they are 



different from the conventional resonant tunneling, because the conventional tunneling is 
completely forbidden for \E\ < A. In fact, these peaks come from the Andreev reflections. 
For the finite-sized pristine nanotubes, the resonant states are symmetric around the Fermi 
energy as mentioned above. Furthermore, the chemical potential of the right superconduct- 
ing lead /i/j = is lined up with the Fermi energy of the nanotube. It is just located in 
the middle of two symmetric states with energy e^ and — e.j. When the electron incident 
from the left lead has the energy e^ corresponding the ith state of the nanotube, a hole can 
propagate back to the state with the energy — Cj. Then a Cooper pair creates in the right 
superconducting lead because of Andreev reflection. The conductance for the nanotube with 
L = 4 is plotted in Fig. 1(c). It is evident that both positions and heights of the resonant 
peaks depend on the nanotube length. There is one resonant peak with the amplitude of 
two units at the Fermi energy, which is absent for the nanotube with L = 3. Similarly, there 
is a resonant Andreev peak at the Fermi energy for the N/CNT/S system with L = 4 as 
shown in Fig. 1(d). This is attributed to the electronic properties of the nanotubes. The 
band structure of armchair nanotubes consists of two non-degenerate bands that cross the 
Fermi level at kp = 27r/3a, with lattice constant a. Finite size effects in carbon nanotubes 
lead to the quantization of the energy levels. In general, one resonant peak appears at the 
Fermi energy with L = 3N + 1 (A^ denoting the number of carbon repeat units), because 
kp = 27r/3a is now an allowed wave vector, a large conductance exists due to a crossing 
of two resonant states at the Fermi energy.— For other lengths, kp is not an allowed wave 
vector and no resonant state exists at the Fermi level, thus conductance is much smaller due 
to the energy gap between the resonant states. These are referred to as on-resonance and 
off-resonance of the Andreev reflection, respectively. 

The impurity, such as the substitutional boron (nitrogen) and the vacancy, can greatly 
change the electronic structure of the nanotubes and then the transport properties. In 
general, the impurity increases the normal reflection.i^ The substitutional boron or nitrogen 
impurity in the infinite carbon nanotube lead to a quasibound state near the lower or 
upper subbands.— However, it is quite different for finite-sized carbon nanotubes. New 
resonant state appears when the incoming electron energies match those of the quasibound 
states induced by the impurity. The substitutional boron effects on the conductance of the 
N/CNT/N system with L = 3 are clearly shown in Fig. 2(a). Compared with Fig. 1(a), 
the impurity leads to one new peak below the Fermi energy. The resonant states for the 



conduction and valence bands are quite different due to the impurity, because the electron- 
hole symmetry, present in perfect nanotubes within the vr-band approximation, is broken.— 
The positions of the resonant peaks is not symmetric around the Fermi energy Ep = due 
to the symmetry broken. The corresponding Andreev reflection probability is plotted in Fig. 
2(b). It can be clearly seen the substitutional boron greatly reduces the Andreev reflection 
probability. Because of the broken of the electron-hole symmetry mentioned above, the 
energy levels of the finite-sized nanotubes are not symmetric around Ep = 0. Because the 
energy level e^ does not have a corresponding one with — e^, the condition for the Andreev 
reflection has been broken. When the electron incident from the left lead has the energy e^ 
corresponding to the ith state, a hole can not propagate back to the state with the energy 
that is different from — e^. The Andreev reflection is then reduced distinctly by the impurity. 
Fig. 2(c) shows the conductance for the nanotube with L = 4. It is seen that the original 
peak at the Fermi energy with the height of two units are splitted into two ones with the 
height of one unit each. One of the two peaks is still at the Fermi energy and the other one 
is above the Fermi level. The reason is that the single impurity breaks the mirror symmetry 
planes containing the tube axis, and then the two resonant states at the Fermi energy are 
now splitted into two ones. Fig. 2(d) shows the corresponding Andreev reflection probability 
for L = 4. The Andreev reflection for L = 4 is stronger than that for L = 3, because the 
resonant states are more symmetric for L = 4. Furthermore, there is also one peak at the 
Fermi energy due to the reasons mentioned above. 

The conductance and the corresponding Andreev reflection probability for the nanotube 
with the substitutional nitrogen are plotted in Fig. 3. A substitutional nitrogen has similar 
effects on the conductance of the N/CNT/N system as the boron, but it induces one resonant 
state above the Fermi energy. The resonant state associated with boron or nitrogen is 
analogous to the acceptor or donor state in semiconductors.— The nitrogen impurity also 
reduces the Andreev reflection, because it breaks the symmetry of the resonant states of the 
finite-sized nanotubes. Because the electron-hole spectra becomes more symmetrical with 
stronger U, the degree of the symmetry broken is slighter. The reduction is not so heavily 
as that caused by the boron. Then the Andreev reflection becomes larger again for stronger 
U. Fig. 4 shows the effects of the vacancy on the conductance and the Andreev reflection. 
For the very large U, as shown in Fig. 4(a), the resonant states become symmetric again. 
However, it is quite different from Fig. 1(a). One new peak emerges near the Fermi energy, 



which is quite different from that of the perfect nanotubes. The position of the resonant 
state caused by the impurity approaches to the Fermi energy for a very strong L'^.— Then, 
it is expected that the vacancy dose not reduce the Andreev reflection, which can be seen 
from Fig. 4(b). For the off-resonant nanotube, the vacancy causes the appearance of a 
resonance state at the Fermi energy, where it is originally zero for the perfect nanotube. For 
the on-resonant nanotube, the vacancy induces a new resonant peak at the Fermi energy as 
shown in Fig. 4(c). The original one at the Fermi energy is splitted into two smaller ones 
around the Fermi energy. Because the electron-hole symmetry is recovered, the Andreev 
reflection becomes stronger again as shown in Fig. 4(d). 

The Andreev reflection current versus the gate voltage is also investigated at different 
bias. As shown in Fig. 5(a)-(c), the current for the nanotube with L = 3 exhibits a single 
series of peaks with different spacings at very small bias, which reflects the band structure 
of the finite-sized nanotubes with the substitutional impurity. At a larger bias, as shown in 
Fig. 5(d)-(e), a series of extra peaks emerges. And they are located in the middle of the 
original ones at the low bias. Furthermore, the amplitude of some original peaks become 
larger and reaches a height three times its value at the low bias. The reason is that more 
energy levels contribute to the Andreev reflection at high bias. When V > Ae/2, a series 
of extra peaks can emerge between the original ones. And when V > Ae, the height of the 
original peaks can become about three times the original value.— Due to the different energy 
spacings, for some high bias, both of the conditions mentioned above can be satisfied. The 
impurity effects on the Andreev reflection current for the tube with L = 4 are clearly shown 
in Fig. 6. There are more peaks in the current at the same bias compared with Fig. 5. For 
nanotube with longer length, the spacing between energy levels become smaller. Because 
more resonant states contribute to the current at the same bias, there are more peaks in 
the Andreev reflection current. The complicated resonance patterns with different kinds of 
peaks depend on the bias, the level spacing of the CNT. It is also found that the exact 
location of the impurity has less influence on this significant feature than the strength of the 
impurity. 

In summary, the probability and the current of the Andreev reflection for the N/CNT/S 
hybrid system are studied in detail. The Andreev reflection exhibits the on-resonance and 
off-resonance behavior at the Fermi energy for the nanotubes with L = 3N + 1 and other 
lengths, respectively. The substitutional boron, nitrogen and vacancy are simulated by 



different impurity strength U. The substitutional boron or nitrogen breaks the electron-hole 
symmetry of the nanotube and reduces the Andreev reflection greatly. The symmetry broken 
and the reduction depend on the impurity strength distinctly. The vacancy with very strong 
impurity strength does not reduce Andreev reflection, because it keeps the electron-hole 
symmetry. It can change the off-resonance to the on-resonance of the Andreev reflection 
at the Fermi energy. The Andreev reflection current shows a complicated behaviour at 
different bias, which reflect the complex band structure of the finite-sized nanotubes with 
the impurity. 
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Figure Captions 

Fig. 1. (a) and (b) are the conductance G and Andreev reflection probability Ta as a 
function of electron energy E for the N/CNT/N and N/CNT/S systems of L = 3, respec- 
tively, (c) and (d) are the corresponding G and Ta for the systems of L = 4, respectively. 

Fig. 2. (a) and (b) are the conductance G and Andreev reflection probability Ta as a 
function of electron energy E for the N/CNT/N and N/CNT/S systems of L = 3 with a 
boron impurity, respectively, (c) and (d) are the corresponding G and Ta for the systems of 
L = 4, respectively. 

Fig. 3. (a) and (b) are the conductance G and Andreev reflection probability Ta as a 
function of electron energy E for the N/CNT/N and N/CNT/S systems of L = 3 with a 
nitrogen impurity, respectively, (c) and (d) are the corresponding G and Ta for the systems 
of L = 4, respectively. 

Fig. 4. (a) and (b) are the conductance G and Andreev reflection probability Ta as a 
function of electron energy E for the N/CNT/N and N/CNT/S systems of L = 3 with a 
vacancy, respectively, (c) and (d) are the corresponding G and Ta for the systems of L = 4, 
respectively. 

Fig. 5. (a)-(c) are the Andreev reflection current Ia as a function of the gate voltage Vg 
at the bias V = 0.03 with a boron, a nitrogen, and a vacancy, respectively, (d)-(e) are the 
corresponding Ia at the bias V = 0.3. Here L = 3. 

Fig. 6. (a)-(c) are the Andreev reflection current I a as a function of the gate voltage Vg 
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at the bias V = 0.03 with a boron, a nitrogen, and a vacancy, respectively, (d)-(e) are the 
corresponding Ia at the bias V = 0.3. Here L = 4. 
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